Source of Acquisition 
NASA Marshall Space Flight Center 


Scaling Equations for Ballistic Modeling of Solid Rocket 

Motor Case Breach 


Joshua E. McMillin* 

ATK Launch Systems, Brigham City, UT, 84302-0707 


This paper explores the development of a series of scaling equations that can take a 
known nominal motor performance and scale it for small and growing case failures. This 
model was developed for the Malfunction-Turn Study as part of Return to Flight activities 
for the Space Shuttle program. To verify the model, data from the Challenger accident (STS- 
51L) were used. The model is able to predict the motor performance beyond the last 
recorded Challenger data and show how the failed right hand booster would have performed 
if the vehicle had remained intact. 
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Nomenclature 

Ratio of the specific heats for a gas 

Propellant density 

Combustion gas density 

45 th Space Wing - Patrick Air Force Base 

Propellant bum rate constant 

Nozzle throat area 

Area of case breach 

Nozzle exit area 

Total propellant surface area 

Total effective throat area 

Effective exhaust velocity 

C* at reference pressure 

Thrust coefficient 

Total bum distance 

Thrust 

Anomalous thrust due to case breach 

Data index used when scaling nominal performance 

Mach number or total expended mass 

Mass flow rate 

Bum rate pressure exponent 

Pressure 

Ambient pressure 
Chamber stagnation pressure 
Nozzle exit plane static pressure 
Reference pressure 
C* sensitivity to pressure 
Propellant bum rate 
Reusable Solid Rocket Motor 
Space Transportation System 
Time 

Time interval 

Value of parameter A for nominal motor 
Value of parameter A for deviating motor 
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I. Introduction 

T HIS paper explores the development of a series of scaling equations that can scale known nominal motor 
performance for small and growing case failures that result in overboard combustion gas venting but not 
immediate case rupture. Previous modeling efforts assumed a simple 2-12 second decay to zero thrust from the time 
of failure. In June of 2003, the RSRM ballistics team was asked to revisit these analyses to support the updated 
Launch Area Risk Assessment performed by the 45th Space Wing (45SW) and ACTA, Inc. It was determined that a 
more realistic, or physics based, failure model was needed. The resulting failure model simulated a growing hole in 
the case wall. The model parametrically scales motor performance for a given time of failure initiation. The failure 
affects time, chamber pressure, thrust, and mass flow rates. This model was provided to NASA’s Johnson Space 
Center and incorporated into the 1 st Stage Malfunction-Turn Study. 

The evaluation of these failure models was purposefully initiated for conducting contingency analyses. They do 
not reflect a true shortcoming of the motor design. Failures such as these are extremely unlikely due to the 
robustness of the motor case and joint seal design. 

II. Background 

The 1st Stage Malfunction-Turn Study provided nominal and malfunction turn trajectories in support of the 
updated Launch Area Risk Assessment performed by the 45SW and ACTA, Inc. The trajectories were used as inputs 
for the calculation of casualty expectation values, which the range safety community used to evaluate public risk. 
Prior to this study, malfunction-turn data had not been provided to the 45SW since STS-5. The delivery of the 
updated malfunction-turn trajectories was coordinated by the Space Shuttle Range Safety Panel. 

The previous (STS-5) Malfunction-Turn Study incorporated two types of RSRM failures modes; single actuator 
failures and thrust decay failures. The thrust decay methodology employed in that study involved linear or 
exponential thrust-decay profiles where the thrust decayed to zero over a 2-12 second time interval. Due to the 
unrealistic nature of these previous models, NASA requested that models be developed that are both more realistic 
and physics based. The new thrust decay models developed by ATK Launch Systems served as inputs to the 
trajectory simulation used to generate the data delivered to the 45 SW. 

The original thrust decay model was intended to simulate a failure where the nozzle throat is burned through. 
Bum through and subsequent ejection of the nozzle are considered unreasonable for the RSRM because the nozzle 
throat cannot be physically ejected in one piece. For this situation to occur, either the throat insert would have to 
break apart, or the entire nozzle would need to be destroyed. The safety factors in the throat region are sufficiently 
high to keep this from occurring. In addition, it is believed that a complete ejection of the nozzle would not 
extinguish the motor. Therefore, the propellant grain would continue to bum and thrust would not decay to zero. 

In an effort to re-evaluate the models used, ATK Launch Systems developed two models for the new study, 
neither failure scenario is considered likely. The first model simulated a failure of nozzle joint #1. The second model 
simulated a growing breach of the pressure chamber. 

Although such a failure is highly unlikely, the failure of nozzle joint #1 was selected as the most realistic and 
readily modeled failure point of the RSRM nozzle. Because it has less material, the exit cone liner has a lower safety 
factor on char and erosion than does the throat, making the potential for bum through at least theoretically higher. 
Nozzle joint #1 is the location where the aft exit cone attaches to the forward exit cone of the nozzle assembly. 
Figure 1 shows a cut away schematic of the RSRM nozzle, with the joint locations identified. Ideally, chamber 
conditions forward of the throat plane are independent of conditions aft of the throat plane, so a failure of joint #1 
would not impact the propellant grain, nor would the motor performance decay to zero thrust. This was modeled 
simply as an instantaneous change in the expansion ratio of the nozzle. The math for this can be found in most 
textbooks covering compressible flow through a converging diverging nozzle. 

The degradation of thrust resulting from the nozzle joint #1 failure resulted in a significantly less severe 
malfunction-turn than the original study. Based on the assumptions of the Malfunction-Turn Study, and barring any 
secondary failures, the STS would still reach staging with the nozzle joint #1 failure. Therefore, it was requested that 
ATK Launch Systems also provide a model that would lead to continuous, more detrimental, thrust decay. For this 
level of thrust decay to occur, the nozzle throat, nozzle assembly, or the pressure chamber must fail in some way. 
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Figure 1: RSRM Nozzle Joint Locations 


Failure of the throat through rapid erosion was considered extremely implausible because of the large safety 
factors and historical performance of the throat ring. Failure of the pressure chamber may occur if a joint seal or the 
pressure vessel was breached in some other way. The case breach scenario is also highly unlikely; but was chosen 
because a joint failure actually happened on STS-51L. The failure model is derived such that it is not limited to the 
joint areas, and could represent a venting breach anywhere in the pressure vessel. Additionally, this same 
methodology could be used, with some alteration, to model high levels of throat erosion. 

Since the Challenger accident, the field joints have been redesigned to prevent an STS-51L type of failure and 
joint heaters have been added to keep the sealing O-rings at temperatures warm enough to ensure the O-rings track 
joint movements with large safety margins. 


III. Parametric Model of Case Breach 


A. Scaling Equations 

Modeling the case breach failure is not as straight forward as modeling a failure of nozzle joint #1. The case 
breach model is intended to simulate an event in which the pressure chamber is compromised, leading to a thrust 
decay scenario. Timing, pressure, mass flow, and thrust are all affected due to the nature of this failure. The case 
breach is modeled as a second nozzle with increasing throat area. The model does not account for axial pressure 
drop, but rather utilizes a 0-D scaling approach to simulate the failure. This assumption allows the model to be used 
regardless of where the case breach is located. The equations scale a nominal performance curve for a given failure 
onset time. Per requirements, the model is able to initiate the failure at any given time. Section B will detail the 
derivation of the scaling equations for time, pressure, thrust, and mass flow. 

Currently the RSRM program uses similar scaling equations for flight predictions. These standard scaling 
equations adjust a nominal “block motor” performance based on loaded propellant weight, predicted bum rate, and 
propellant mean bulk temperature. The scaling equations derived for the case breach failure can be used in 
conjunction with these standard scaling equations. 

Figures 2, 3, and 4 show the results of this model applied to a nominal RSRM with a failure initiation time of 25 
seconds. The failure is modeled as a circular area with a constant radial growth rate. These results are given as an 
example only. Final use of the model by engineers at NASA Johnson Space Center used multiple failure initiation 
times and locations. 
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Figure 3: Results From the Parametric Model for Thrust 


Malfunction Analysis - Mass Flow Rates 
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Figure 4: Results from the Parametric model for Mass Flow 


B. Deriving the Scaling Equations 

A complete list of nomenclature 
used in the derivation of these 
equations can be found at the 
beginning of this paper. 

The key assumption in the 
scaling equation approach is that 
each data point in the nominal data 
file represents a fixed bum distance 
or a fixed propellant geometry. 
Although the breach is growing, the 
analysis assumes that steady state 
equations are still applicable. The 
result is a set of data correlating time, 
chamber pressure, thrust for both the 
nozzle and the breach, mass flow for 
both the nozzle and the breach, and 
the breach area. 

1. Pressure Scaling 

Using steady state equations, the 
nominal head-end pressure can be 
scaled as a function of the nominal 
and deviant effective throat areas to 
predict the deviant pressure trace. 
The nominal throat area (A * nom ) is the 
original area of the nozzle throat 
associated with a given data point. 
The deviant throat area (A l dev ) is the 
combined total area of the new 
nozzle throat and case breach. The 
deviant nozzle throat may be slightly 
different from the nominal nozzle 
throat because of normal throat 
erosion occurring over a longer 
period. 

The steady state assumption 
means that mass flow being 
generated at the burning propellant 
surface is equal to the mass flow 
being expelled through the nozzle 
and case breach. There is no 
appreciable change in the mass 
stored in the chamber volume. 
Assuming each data point in the 
nominal performance represents a 
fixed bum distance, and therefore a 
fixed surface area ( A s ), the surface 
area term cancels out along with the 
propellant density (p pr op) and bum 
rate constant term (a). 

The pressure scaling equation 
begins with the fundamental 
equations for Bum rate, Eq. (1), and 
mass flow, Eq. (2). Ultimately, the 
result is Eq. (8). 
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Equation (1) is a general formula for estimating the bum rate of a solid propellant as a function of pressure (P) 
and empirically derived constants a and n. 


r = aP n 


( 1 ) 


Equation (2) calculates the mass flow from the propellant surface as a function of surface area the density of the 
propellant and the bum rate. 


™ = A sP proof 


( 2 ) 


Equation (3) is the definition for the Characteristic Exhaust Velocity (C*), which relates pressure, total throat 
area (A% and mass flow. 

C = PA' / m (3) 

Equation (4) is the definition for the C* Sensitivity to Pressure ( q ), an empirically derived constant. 

q = d ln(c* y</ ln(p) (4) 

The equation for C* can be rewritten using Eqs. (1) through (3) to be a function of pressure and surface area. 

c = P'-A'/(A s p r „„a) (5) 

Equation (6) uses Eq. (4) to relate C* at different pressure levels. Equation (7), which comes from substituting 
Eq. (5) in to Eq. (6), can be solved to give the pressure for the deviant motor, Eq. (8). 


C\r 4 = c-Jp^ipJ - c;Jp^ f /p,J 


( 6 ) 
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(7) 


( 8 ) 


2. Time Scaling 

To create scaling equations that are simple and easy to solve, an approximation was made. This approximation 
was that the pressure, nozzle throat area, and breach area remain sufficiently constant between two adjacent points in 
the nominal performance data. The longer the failure was allowed to grow the more error this assumption 
introduced. This error was insignificant for sufficiently small increments of time. As seen in section 5 this 
assumption also affected the mass flow calculations. A more accurate derivation is presented in section 6 

The key assumption for the scaling equations is that each data point in the nominal performance data represents a 
fixed bum distance (£>), or a fixed propellant geometry. Therefore, the distance burned between points [M] and [/] 
remains constant as in Eq. (10). 


D = r At 

(9) 

D =D, 

nom dev 

(10) 
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Combining Eq. (9), Eq. (10), and Eq. (1) relates the nominal pressure and time delta (A t nom ) with the predicted 
deviant pressure and time delta (A t dev ) as seen in Eq. (11). 


aP n At = aP" At, 

nom nom dev dev 


( 11 ) 


Substituting Eq. (8) into Eq. (1 1), the deviant time delta between two data points can be expressed by Eq. (12) 


^ dev ^ nom T nom / 4fcv ) 


\-n!(\-n-q) 


( 12 ) 


The time for the deviant performance (T devi ) is the deviant time delta plus the deviant time value for the previous 
data point ( T deVti .;). 


^devd ^devj - 1 ^ dev 


(13) 


The area of the breach (A B ) can be modeled as any function of time. For the purpose of the Malfunction-Turn 
Study, this was modeled as a simple circular hole with a fixed growth rate beginning at a given time (7^). 


A,=f(T d „-T,) 


(14) 


The total effective throat area (A r dev ) is the sum of the breach area (A B ) and the nozzle throat area (A* rdev ). 


A dev ~ A r dev + A b 


(15) 


Because the breach area increases with increasing time, and because of the slow erosion of the nozzle throat, it is 
easiest to iterate on Eqs. (12) - (15) until a convergence on time is reached. The spreadsheet provided by ATK 
Launch Systems to engineers at NASA Johnson Space Center for the Malfunction-Turn Study converges to ±0.001 
seconds at all times within eight iterations for failure onset at initiation. The bigger the failure is allowed to grow, 
the more iteration is needed to converge within a given tolerance. 

3. Nozzle Thrust Scaling 

The nozzle thrust (F Nozz ie) scaling is obtained by assuming that for a given data point in the nominal performance, 
the thrust coefficient, defined in Eq. (16), remains constant as in Eq. (17). In effect, this also states that the 
expansion ratio is unaffected. The small change in nozzle expansion ratio is neglected because of the added 
modeling complexity and because of its relatively small effect. 


r = F 

^ F Nozzle 


/ PA * 


(16) 


c =c 

^ F,dev ^ F,nom 


(17) 


Substituting Eq. (16) into Eq. (17) results in Eq.(18). This relates thrust to the product of pressure and throat area 
for both the nominal and deviating conditions. 


F Nozzle, dev j i/^dev ^T dev ) ^ Nozzle, nom / i/^nom ^T nom ) 


(18) 


Equation (19) scales thrust from the nominal performance as a function of effective throat areas using Eq. (8) 
with Eq. (18). 


F = F 

Nozzle, dev Nozzle, nom 


( A r Jev ! A T„ om \ A T„ om / A dev) 


i/(\-q-n) 


(19) 
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4. Breach Thrust 

The derivation for the breach thrust (F B ) begins with the ideal thrust equation 1 , Eq. (20). This can be greatly 
simplified given the assumptions of vacuum ambient conditions and choked flow through an orifice, where the 
throat (A*) and exit ( A e ) areas are both equal to the breach area (A B ). 


F - ^•/>«r[{2/0'-l)}{2/(y +1)} 0 - +(P, -P.)A 

a,=a-.a, 

Equation (22) is the isentropic relationship between Mach number and pressure. 

PjP c = ^ + 0.5{y-\)M 2 Y (M) 


( 20 ) 

( 21 ) 

( 22 ) 


The Mach number for choked flow through an orifice is 1.0 and Eqs. (20) - (22) combine to yield Eq. (23). 

F t = p ^ a b (r + iX2/(r + i)y /0 " 1) < 23) 

Adjustments for non-ideal flow losses or dimensional effects are not considered. For the purposes of the 
Malfunction-Turn Study, the off nominal thrust generated by the breech is kept as high as possible to maximize 
thrust imbalance between the failed and nominal boosters. 

5. Mass Flow Rate Scaling 

Mass flow rate scaling is based on conservation of mass, Eq. (24), and assumes that the difference in mass 
storage is negligible. The mass flowed from data point to data point must be the same in order to meet the 
assumption that each data point represents a fixed bum distance and fixed propellant geometry. Using a numerical 
integration technique (trapezoidal rule), Eq. (25), the total mass flow is determined using Eq. (26). 


^ i,dev ^ i,nom 


M t = M m + 0.5(m. + m M )A t 


(24) 

(25) 


Equations (24) and (25) combine to form Eq. (26). This is the equation used in the Malfunction-Turn Study to 
scale the nominal mass flow data for a malfunctioning booster. 


=(<„„« + ) - ^,-Uev 


(26) 


Equations (29) and (30) are used to determine the mass flow being ejected from the nozzle and case breach 
respectively. These were based on the assumption that the same gas species, at identical conditions (constant gas 
density, p g ), are flowing through both the nozzle throat and the breach. Axial pressure drop is ignored in the scaling 
equations to simplify the approach and allow simulations to use the same equations for a failure at any location. 
With this assumption, the individual mass flow rates become proportional to the breach and nozzle throat areas. 


m/A ‘ 


m = p g V ■ A 

(27) 

— ^Nozzle! A = Wl B / A b 

(28) 

™ Nozzle =rh i A'/A' 

(29) 

=M i A B /A' 

(30) 
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6. Corrected Mass Flow and Time Scaling 

From Eqs. (1) (2) and (8) listed previously, it is also possible to calculate Eq. (31). Unfortunately, with the time 
scaling derived in section 2, this equation will violate the conservation of mass requirement. On an instantaneous 
basis, this difference is small; however, integrated over time the error became evident. For failure onset at initiation 
(T=0) the error accumulated to almost 0. 15% of the RSRM total expended mass. 

<e w = < ld X ( A i / A L • ( 31 ) 

Upon further investigation, it was found that this error became smaller and smaller with decreasing the 
magnitude of time increments used in the nominal performance data. Ideally, Eqs. (26) and (31) should yield the 
same results. The difference has been identified as an inadequacy in Eq. (9), which assumes that the pressure change 
between data points in the nominal performance is negligible. This assumption breaks down as the failure progresses 
and the change in the time steps and pressure become larger. It also breaks down for motors having more radical 
changes in pressure. Equation (32) is a more correct equation for the distance burned between two points. 

D = \rdt (32) 

Equation (12) can be replaced using Eq. (1), Eq. (10), and a trapezoidal integration rule to approximate Eq. (32) 
°- 5fl fev + P'-Uev Kto, = °M P , n no m + Pi-X.nom Kom (33) 


This can be written in terms of the effective throat and known pressures. The new equation is somewhat longer 
than Eq. (12), but it is easy enough to translate into a computer program or spreadsheet. 


At. = At 

dev nom 


i A t w(i-»-?) (p ip yb 1 

f/ 1 i.nom / ^i,dev ) ^ V i-\,dev / r i,nom J J ~ t ' 

I(a* I a 1 Y^~ n ~ q Hp. Ip. V + (p Ip V \ 1 

_ i,dev J V i,nom / 1 i-\,nom ) ' Vr i-\,dev / 1 i-\,nom ) j 


(34) 


Substituting Eq. (34) for Eq. (12) and Eq. (31) for Eq. (30) results in a significantly reduced amount of error in 
the total integrated mass flow (<0.0003%). 

C. Modeling the Case Breach 

For the purpose of the malfunction study, the case breach was modeled as a circular opening in the pressure 
vessel with a fixed growth rate. There are no known design deficiencies, or preventable scenarios, that could lead to 
this failure. Modeling of the case breach was performed to provide information used in contingency planning for 
unforeseeable events. The shape of the failure was selected based on ease of modeling. The growth rate was based 
on a charring material ablation model of hot propellant combustion products flowing through a circular hole in a flat 
steel plate. This was considered conservative because it did not include any protective materials such as insulation, 
liner, or propellant. STS-51L data was also evaluated to confirm the appropriateness of the results. Analysis of the 
available pressure data indicates that the STS-51L failure grew at a slower rate. Therefore, the model used in the 
Malfunction-Turn Study was conservative and generated a larger thrust imbalance. 

D. STS-51L 

The STS-51L case breach was initiated by a failure of the aft field joint of the right hand booster. The case 
breach failure model was used to simulate the performance of that right hand motor. A time of 58.788 seconds was 
picked as the failure initiation time because it was the time determined by the Challenger accident investigation 2 for 
first evidence of flame. The joint failure was modeled in two different ways. First as a circular hole growing with 
constant radial growth rate. A radial erosion rate of 0.25-in/sec (Fig. 5) best matched available pressure data with 
this assumption. A radial erosion rate of 0.32-in/sec (Fig. 6) was needed to match the final failure area with a 
constant radial growth. Second, the failure was modeled with a constant area growth rate. With this assumption, a 
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STS-51L Model Comparison 



Figure 5: STS-51L Modeled with 0.25-in/sec Radial Growth Rate 


STS-51L Model Comparison 


| LH Booster — RH Booster Breach Model — HPM Nominal* | 



Figure 6: STS-51L Modeled With 0.32-in/sec Radial Growth Rate 


STS-51L Model Comparison 


j LH Booster RH Booster — Breach Model HPM Nominal* | 



constant area growth rate of 3-in 2 /sec best 
matched the available chamber pressure 
data. Clearly, the results obtained depend 
on the assumptions made about the 
geometry and growth of the venting 
failure. 

The recovery of the STS-51L 
hardware allowed for an estimate of the 
final dimensions of the case failure. A 
bottom portion of the aft cylinder of the 
aft center segment was recovered with a 
burned out area approximately 15 axial 
inches and 28.25 circumferential inches. 
A portion of the aft attach cylinder was 
also recovered with evidence of burn 
through defining the bottom edge of the 
STS-51L failure. Figure 8 and Fig. 9 show 
the dimensions of the case segments 
recovered with evidence of the bum 
through. From these two pieces, it was 
possible to estimate the size of the failure 
at the time of destruction. According to 
the failure investigation 3 , the failure was 
approximately 862-in 2 . If the failure can 
be modeled as a circular failure with 
constant radial growth rate, the resulting 
rate based on an area of 862-in 2 and 
destruction time of 110.250-seconds 
(51.462-seconds after first flame) 2 would 
be approximately 0.32-in/sec. Results of 
the parametric model using this erosion 
rate are shown in Fig. 6. The Circular 
assumption does not initially decay 
pressure as quickly as the measured data 
indicates. Assuming an erosion rate of 
0.32-in/sec causes the predicted pressure 
in excess of the measured data after less 
than ten seconds. 

To best match the measured data, a 
constant area growth rate of 3-in 2 /sec, as 
shown in Fig. 7, was determined. This 
growth rate does not lead to the size of the 
apparent failure at the time of motor 
destruction. To achieve the observed 
failure size in the time from first flame to 
destruction, with a constant area growth 
rate would require a rate of 16.75-in 2 /sec. 
These differences highlight the fact that 
the failure itself was not so easily 
modeled. Reality proves that the growth 
rate of the failure may not be linear, or 
even continuous. In addition, events after 
the loss of the data stream may have 
caused additional damage. 


Figure 7 STS-51L Modeled with failure growth rate of 3-in 2 /sec 


9 

American Institute of Aeronautics and Astronautics 






Ml . IMHII 


Figure A 8 Aft Center Segment Right Hend SRB Contact No. 131 

Figure 8: Contact 131, RH Aft Center Segment Aft Cylinder 4 



Figure A -22 All Sogmam Right Hand SRB Contact 712 

Figure 9: Contact 712, Aft Segment Right Hand Booster Forward Cylinder 4 
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IV. CONCLUSIONS 

The scaling approach can be used to model motor performance for a given failure. Using fundamental 
engineering equations, a robust algorithm has been derived for modeling the performance of a booster with a breach 
in the pressure vessel. For future analyses, simple modifications can be made to this model. Such changes may 
include modifying the failure growth rate, failure shape, or using the scaling approach to model high rate nozzle 
throat erosion that would also require modeling the change in expansion ratio. Evaluation of the Challenger data 
shows that the model is reasonable and adaptable. This tool provides the ability to model abnormal motor 
performance for use in evaluating contingency scenarios. 

REFERENCES 

1 Humble, R. W., Gary N. H., and Wiley J. L., Space Propulsion Analysis and Design , 1 st ed.-Revised, The McGraw-Hill 
Companies, Inc., New York, 1995, Chap. 3.3. 

2 “Report of the PRESIDENTIAL COMMISSION on the Space Shuttle Challenger Accident,” Vol. 1, Chap. 3, Washington 
D.C., 1986. 

3 “Report of the PRESIDENTIAL COMMISSION on the Space Shuttle Challenger Accident,” Vol. 3, Appendix O, 

Enclosure 9, Washington D.C., 1986. 

4 “Report of the PRESIDENTIAL COMMISSION on the Space Shuttle Challenger Accident,” Vol. 2, Appendix L, Appendix 
A, Washington D.C., 1986. 


11 

American Institute of Aeronautics and Astronautics 



